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ABSTRACT 

We use a semi-analytic model for disk galaxies to explore the origin of the time evo- 
lution and small scatter of the galaxy SFR sequence — the tight correlation between 
star-formation rate (SFR) and stellar mass (Mgtar)- The steep decline of SFR from 
z ~ 2 to the present, at fixed Mstar, is a consequence of the following: First, disk galax- 
ies are in a steady state with the SFR following the net (i.e., inflow minus outflow) 
gas accretion rate. The evolution of the SFR sequence is determined by evolution in 
the cosmological specific accretion rates, oc (1 -I- z)^-^^, but is found to be independent 
of feedback. Although feedback determines the outfiow rates, it shifts galaxies along 
the SFR sequence, leaving its zero point invariant. Second, the conversion of accretion 
rate to SFR is materialized through gas density, not gas mass. Although the model 
SFR is an increasing function of both gas mass fraction and gas density, only the gas 
densities are predicted to evolve significantly with redshift. Third, star formation is 
fueled by molecular gas. Since the molecular gas fraction increases monotonically with 
increasing gas density, the model predicts strong evolution in the molecular gas frac- 
tions, increasing by an order of magnitude from z = to z ~ 2. On the other hand, 
the model predicts that the effective surface density of atomic gas is ^ IOMqpc"^, 
independent of redshift, stellar mass or feedback. Our model suggests that the scatter 
in the SFR sequence reflects variations in the gas accretion history, and thus is insen- 
sitive to stellar mass, redshift or feedback. The large scatter in halo spin contributes 
negligibly, because it scatters galaxies along the SFR sequence. An observational con- 
sequence of this is that the scatter in the SFR sequence is independent of the size 
(both stellar and gaseous) of galaxy disks. 

Key words: galaxies: evolution - galaxies: formation - galaxies: fundamental pa- 
rameters - galaxies: haloes - galaxies: high redshift - galaxies: spiral 



1 INTRODUCTION 

Understanding the star formation history of the universe is a 
major goal of modern astronomy. Starting at redshift z = Q, 
the comoving cosmic star formation rate (SFR) density (in 
units of Mq yr~^ Mpc~^) rises by an order of magnitude to 
z ~ 1 (e.g. Lilly et al. 1996; Schiminovich et al. 2005; Hop- 
kins & Beacom 2006; Villar et al. 2008, Sobral et al. 2009), 
reaches a peak around z ~ 2 and declines to higher redshifts 
(Madau et al. 1998; Perez- Gonzalez et al. 2008; Bouwens 
et al. 2008). 
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A number of studies have shown that the decline in the 
comoving cosmic SFR density from 2 ~ 1 to the present 
is not due to a decline in the major merger rate (Bell 
et al. 2005; Wolf et al. 2005; Melbourne et al. 2005; Lotz 
et al. 2008; Sobral et al. 2009). Rather, since the majority 
of all stars seem to form in galactic disks (at least below 
2: ~ 2), it is generally believed that the drop in comoving 
cosmic SFR density from z ~ 1 to the present is related 
to 'quiescent' processes, e.g., gas consumption and/or a re- 
duced rate of accretion of fresh gas from the cosmic web. 

Great progress is currently being made in measuring 
molecular gas masses in high redshift galaxies (e.g. Tacconi 
et al. 2006, 2008, 2010; Daddi et al. 2008, 2009b). These 
studies indicate that the ratio between molecular gas and 



© 2009 RAS 



2 Button et al. 



1 — I — I — r 



T 



2 - 



1 - 



/ 0' 




scotter 



9 10 11 

IOg,0 ^s^or [Mo] 



0.5 - 



^ 
cr. 

Li_ 
CO 

o -0.5 

o 



1 - 



M 



star 



:1010 



t 



I I 

;i I I 



-A 



! ! ! -I 

J_ J_ J_ L 

Gonzolez et ol. 2009 
Stork et ol. 2009 
Doddi et o. 2009 
Doddi et cl 2007 
Elboz et ol. 2007 
Noeske et ol. 2007 
Solim et cl. 2007 
Elboz et ol. 2007 
. I . . . I . . . 



0.2 0.4 0.6 



0.8 



Figure 1. Summary of observations from redsliifts 2 ~ to 2 ~ 7. Left-liand panel: SFR sequence at different redshifts. Right-hand 
panel: redshift evolution of the specific star formation rate (SSFR=SFR/Mstar) at a stellar mass of Mstar = 10^° Mq. The horizontal 
error bars indicate the redshift range of the data points, while the vertical error bars indicate the estimated intrinsic scatter in SSFR. 



stellar mass was higher in the past. However, uncertainty in 
the conversion factor between CO luminosity and molecular 
gas mass limits the current accuracy of these measurements. 
In addition there are at present no direct observations of 
atomic gas masses or gas accretion rates in high redshift 
galaxies, and thus the origin of the decline in SFRs since 
2 ~ 2 is not yet understood. 

A step forward in the characterization of the evolution 
of galaxy SFRs was the discovery of a relatively tight cor- 
relation between galaxy SFR and stellar mass (Mstar), with 
a slope just below unity and an intrinsic scatter < 0.3 dex 
(Noeske et al. 2007a; Elbaz et al. 2007; Daddi et al. 2007; 
Salim et al. 2007). Throughout this paper we refer to this 
correlation as the SFR sequence. The zero point appears 
roughly constant from 2 ~ 7 to 2; ~ 2 (Daddi et al. 2009a; 
Stark et al. 2009; Gonzalez et al. 2009), and evolves to lower 
SFR from 2 ~ 2 to 2 = (Daddi et al. 2007). See Fig.[T]and 
Table[T]for a summary of the observational data. The scatter 
of the SFR sequence appears to be invariant with redshift, 
while there is some indication that the slope is shallower at 
lower redshifts. 

The two main questions that we address in this paper 

axe: 

(i) What is the physical cause of the decline in the specific 
star formation rate (SSFR=SFR/Mstar) since 2 ~ 2? 

(ii) Why does the SFR sequence have such a small scat- 
ter? 

A plausible explanation for the decline in SSFR since 
2 ~ 2 is a decline in the gas accretion rate onto galaxies. 
For dark matter haloes the specific accretion rate scales as 
(1 + zf-^^ for 2 < 2 (Birnboim, Dekel, & Neistein 2007). 
Furthermore, scatter in the gas accretion rates from cosmo- 
logical simulations is about 0.3 dex (Dekel et al. 2009). Thus 
the decline in SSFR as well as the small scatter could be 



driven entirely by the properties of gas accretion. However, 
for this to be the case, a mechanism is needed to explain 
why the star formation rate should follow the gas accretion 
rate. 

To get some insight into this problem, we assume that 
the star formation rate surface density follows a Schmidt- 
Kennicutt (SK, Schmidt 1959, Kennicutt 1998) law: 

SspR = espEgjjs, (1) 

and that the gas surface density follows an exponential dis- 
tribution: 

Egas(i?) = Eo exp(-i?/7?gas), (2) 

where Eo is the central surface density of the gas, and i?gas 
is the scale length of the gas disk. The SFR is then given by 

/> 00 

SFR = 27r / EsFRiR)RdR (3) 
Jo 

= espMgasn-'Er' (n > 0) (4) 

Thus, for n > 1, the mechanism which converts higher ac- 
cretion rates into higher star formation rates is either an 
increase in the gas masses or an increase in the gas surface 
densities. 

There is a significant scatter (~ 0.14 — 0.2 dex) in the 
sizes of galactic stellar disks at a given stellar mass or lumi- 
nosity (e.g. Shen et al. 2003; Courteau et al. 2007), which 
corresponds to a scatter in disk surface density of ~ 0.28—0.4 
dex. If the surface density of gas is correlated with that in 
stars, then the observed amount of scatter in the surface 
densities of stellar galactic disks would imply a scatter of 
0.14-0.2 dex in SSFR (assuming n — 1.5), and even larger 
scatter for a higher n. Thus galaxy size could be an impor- 
tant contributor to the scatter in the SFR sequence. 

In this paper we use the semi-analytic disk galaxy evo- 
lution model of Dutton & van den Bosch (2009) to study the 
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Table 1. Summary of the observed SFR sequence from low to high redshift. The relation is of the form SFR = SFRq (A/star /1O1''A/0)='°p<=. 
All stellar masses and SFRs have been converted to a Chabrier IMF. Parameters measured by ourselves are indicated by a *. 



Redshift 


Redshift range 


SFRo 


slope 


Mass range 


original IMF 


Reference 






[Moyr-i] 




log 10 [Mq] 






0.06 


0.015-0.10 


1.29 
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evolution and scatter of the SFR sequence. In §2 we give a 
brief overview of the galaxy evolution model. In §3 we dis- 
cuss the evolution of the zero point of the SFR sequence, and 
whether it is driven by gas accretion rate, gas fractions or 
gas density. The scatter in the SFR sequence is addressed in 
§4, while §5 discusses why the SFR sequence is independent 
of feedback. We summarize our results in §6. 



2 DISK GALAXY EVOLUTION MODEL 

Here we give a brief overview of the disk galaxy evolution 
model used in this paper. A more detailed description of this 
model is given in Button & van den Bosch (2009; hereafter 
BB09). 

The key difference with traditional semi-analytic mod- 
els (e.g. Kauffmann et al. 1993; Cole et al. 1994,2000; 
Somerville & Primack 1999; Hatton et al. 2003; Cattaneo 
et al. 2006; Croton et al. 2006) is that in this model the in- 
flow (due to gas cooling), outflow (due to supernova driven 
winds), star formation rates, and metallicities, are all com- 
puted as a function of radius, rather than being treated as 
global parameters. This allows us to follow self consistently 
the growth of gaseous and stellar galactic disks, and to split 
the gas disk into atomic and molecular phases. Our model 
should thus be considered in the context of previous models 
of this type (e.g. Kauffmann 1996; Avila-Reese et al. 1998; 
Firmani & Avila-Reese 2000; van den Bosch 2001,2002). 

The main assumptions that characterize the framework 
of our models are the following: 

(i) Bark matter haloes around disk galaxies grow by the 
smooth accretion of mass which we model with the Wechsler 
et al. (2002) mass accretion history (MAH). The shape of 
this MAH is specified by the concentration of the halo at 
redshift zero; 

(ii) Bark matter haloes are assumed to be spherical and 
to follow a density distribution given by the NFW profile 
(Navarro, Frenk, & White 1997), which is specified by two 
parameters: mass and concentration. However, halo concen- 
tration is known to be correlated with halo mass, and we use 



mass dependence and evolution of the concentration param- 
eter given by the Bullock et al. (2001a) model with param- 
eters (Maccio et al. 2008) for a WMAP 5th year cosmology 
(Bunkley et al. 2009); 

(iii) Gas that enters the halo is shock heated to the virial 
temperature, and acquires the same distribution of specific 
angular momentum as the dark matter (van den Bosch 
et al. 2002). We use the angular momentum distributions 
of the halo as parametrized by Bullock et al. (2001b) and 
Sharma & Steinmetz (2005); 

(iv) Gas cools radiatively, conserving its specific angu- 
lar momentum, and forms a disk in centrifugal equilibrium. 
Cooling rates are computed using the cooling functions of 
Sutherland & Bopita (1993) under the assumption of coUi- 
sional ionization equilibrium; 

(v) Star formation is assumed to take place in dense 
molecular gas. In order to compute the density of molec- 
ular gas, we follow Blitz & Rosolowsky (2006) and relate 
the molecular gas fraction to the midplane pressure in the 
disk. For completeness, a detailed description of our star for- 
mation prescription is given in Appendix A, where we show 
that it yields a SK star formation law that is truncated at 
low Sgas, in good agreement with observations; 

(vi) Supernova (SN) feedback re-heats some of the cold 
gas, ejecting it from the disk and halo. We assume that 
ejected gas is lost forever. We consider two feedback models, 
one driven by energy (e.g. Bekel & Silk 1986), the other by 
momentum (e.g. Murray, Quataert, & Thompson 2005). For 
both models we assume that the outflow moves at the local 
escape velocity of the disk-halo system, Vesc{R)- The key 
difference between these two models is the dependence of the 
mass loading factor, i] (ratio between outflow mass rate and 
star formation rate), on wind velocity. For the energy driven 
wind rj oc V^c^, whereas for the momentum driven wind 77 oc 
V^^ . Note that since the escape velocity typically varies with 
galacto-centric radius, so do does the mass loading factor. 

(vii) Stars eject metals into the inter-stellar medium, en- 
riching the cold gas; 

(viii) Bruzual & Chariot (2003) stellar population syn- 
thesis models are convolved with the star formation histo- 
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ries and metallicities to derive photometric fluxes in different 
pass bands. 

Each model galaxy consists of five mass components: dark 
matter, hot halo gas, disk mass in stars, disk mass in cold 
gcis, and ejected gas. The dark matter and the hot gas are 
assumed to be distributed in spherical shells, the disk mass 
is assumed to be in infinitesimally thin annuli. Throughout 
this paper we refer to R as radius, t as time (where t = is 
defined as the Big Bang) and z as redshift. 

For each galaxy we set up a radial grid with 200 bins 
quadratically sampled from between 0.001 and 1 times the 
redshift zero virial radius. As described in DB09 for the pur- 
pose of conserving angular momentum it is more convenient 
to use a grid in specific angular momentum, j. We convert 
the initial grid in R, to a grid in j using j'^/G — RM{R), 
where M{R) is the halo mass enclosed within a sphere of 
radius R. We follow the formation and evolution of each 
galaxy using 400 redshift steps, quadratically sampled from 
z = 10 to 2; = 0. For each time step we compute the changes 
in the various mass components in each radial bin. 

2.1 Limitations of the Model 

Our model makes several simplifying assumptions about the 
physical processes that occur during disk galaxy formation 
(mergers, gas accretion, star formation, feedback, adiabatic 
contraction). In some cases it is because the correct physics 
is unknown, in other cases it is because modeling all of the 
details adds extra complication to the model without nec- 
essarily improving our understanding. A discussion of these 
assumptions and the limitations they imply to our model are 
given in DB09. Here, and in Appendix B, we focus on the is- 
sue of gas accretion, and specifically the differences between 
the classical "hot mode" vs. the new "cold mode" accretion 
models. 

Our assumption about the way gas is accreted into 
galaxies (by a cooling flow of shock heated has) is likely in- 
correct. Simulations suggest that disk galaxies accrete most 
of their mass through cold flows, and that in the absence 
of extra heating or outflows the baryon fraction of galax- 
ies is close to the universal value (Birnboim & Dekel 2003; 
Keres et al. 2005; Dekel & Birnboim 2006). As we discuss in 
Appendix B, for haloes with masses below ~ IQ^^Mq, the 
cooling time is short compared to the free fall time, and thus 
the accretion rates onto the central galaxy are the same in 
the classical cooling model as in the cold mode model. The 
key difference between these two accretion models occurs for 
high mass Mvir ^ 10^^ Mq haloes at high 2; > 2 redshifts. In 
this regime the accretion rates are low in the classical cool- 
ing model because cooling is inefficient, but the accretion 
rates can be high in the cold mode model because accretion 
rates are dominated by cold streams penetrating the hot 
halo (Dekel et al. 2009). Thus our results for the cold gas 
accretion rates onto galaxies with stellar masses less than 
about 5 X 1O^°M0 are not expected to be effected by us 
ignoring the phenomenon of cold streams. 

It is possible that the angular momentum acquired by 
galaxies may be different in the hot mode and cold mode 
accretion models. However, non-radiative cosmological hy- 
drodynamical simulations (e.g. van den Bosch et al. 2002; 
Sharma & Steinmetz 2005) have shown that the specific an- 



gular momentum of the gas and dark matter halo are very 
similar, which can be understood from the fact that both 
have experienced the same torques. Thus as long as most of 
the baryons accrete onto the central galaxy (as is the case 
for galaxies in haloes below ~ 10^^ Mq) then the available 
angular momentum is expected to be the same in the hot 
mode and cold mode accretion models. 

2.2 Overview of model parameters 

In this paper we consider three models that only differ 
in their treatment of feedback and the amount of angular 
momentum loss. In addition to a model without feedback 
(NFB), we consider a model with energy driven feedback 
(EFB) and one with momentum driven feedback (MFB). 
The parameters of these models have been tuned as little 
as possible. The NFB model has no tuning. The MFB and 
EFB models have their respective feedback efficiencies and 
amount of angular momentum loss tuned to reproduce the 
zero points of the rotation velocity - stellar mass and stel- 
lar disk size - stellar mass relations at redshift zero. For the 
MFB model this results in maximal feedback and no angu- 
lar momentum loss, for the EFB model this results in 25% 
feedback efficiency and 30% angular momentum loss. De- 
tails of these models, as well as the resulting galaxy scaling 
relations (rotation velocity - stellar mass, stellar disk size 

- stellar mass, gas fraction - stellar mass, gas metallicity - 
stellar mass) and disk structural properties, are discussed 
in DB09 and Dutton (2009). As a summary, the MFB and 
EFB models reproduce the scaling relations reasonable well, 
while the NFB model does not reproduce any of the scaling 
relations. 

The input parameters of our models are as follows. 

(1) Cosmology: f2m,o, ^a, fib, crs, h, n. In this paper we 
adopt a flat ACDM cosmology motivated by the 5th year 
WMAP results (Dunkley et al. 2009). In particular, we adopt 
fim,o = 0.258, = 0.742, fib = 0.044, as = 0.80, h = 0.7, 
and n — 1. 

(2) Dark Halo structure: K, F, aiogc- We assume that 
the concentrations of dark matter haloes of a given mass, 
Al, follow a log-normal distribution with median c{AI) from 
the Bullock et al. (2001a) model, parametrized by _F = 0.01 
and K = 3.7, and scatter aiogc = 0.11. Whenever this 
model yields c < K we set c = K. This parametrization 
nicely reproduces the distribution of halo concentrations of 
relaxed dark matter haloes in cosmological N-body simula- 
tions (e.g. Wechsler et al. 2002; Maccio et al. 2007, 2008; 
Zhao et al. 2009). 

(3) Angular momentum distribution: A, uinA, ce, o"ioga- 
We assume that the spin parameters of dark matter haloes 
of a given mass, M, follow a log-normal distribution with 
median A and scatter ai^x = 0.5 (Bullock et al. 2001b). For 
our models with no outflows and momentum driven out- 
flows we adopt the cosmological value for relaxed haloes 
A = 0.035 (Maccio et al. 2008), while for our model with 
energy driven outflows we assume that 30% of the angular 
momentum has been lost, so that A = 0.025. This choice of 
A is made to reproduce the zero point of the stellar disk size 

- stellar mass relation (see DB09), and is motivated by the 
fact that baryons tend to lose specific angular momentum 
in hydrodynamical cosmological simulations (e.g. Navarro 
& Steinmetz 2000; Piontek & Steinmetz 2009). We assume 
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Figure 2. The SFR sequence for our EFB model with energy driven outflows. The model galaxies are represented by grey points, the 
black error bars show the error on the median SFR in A/star bins of 0.2 dex width, and the black lines enclose 68.3% of the model galaxies. 
The blue dashed lines show the observed relations (only plotted over the mass range for which observations exist), with the dotted blue 
lines showing the observed scatter of 0.3 dex. The red lines show a fit to the model galaxies in the range 9 ^ logio Afstar/[Af0] ^ 11.0. 
The parameters of these fits are given in the lower right corner of each panel. The slope and scatter of the model SFR sequence is roughly 
independent of redshift, whereas the zero point evolves strongly. 



that the specific angular momentum distribution of both 
dark matter and gas in a halo of given mass and spin pa- 
rameter is given by the fitting function presented in Sharma 
& Steinmetz (2005), which is parametrized by a single pa- 
rameter, a. Based on their numerical simulations, we as- 
sume that the distribution of a is log-normal, with a me- 
dian a = 0.9 and scatter aiaga ~ 0.11. We note that the 
parametrization of the angular momentum distribution by 
Bullock et al. (2001b), which is specified by fi, results in 
very similar disk density profiles as the Sharma & Stein- 
metz (2005) parametrization, except for the extremes of the 
distribution of a and (see Button 2009). 

(4) Star formation: esp. Our star formation prescription 
(see Appendix A), is parametrized by a single free param- 
eter, esF, which sets the normalization of Espr at a given 
surface density of dense molecular gas. We tune this param- 
eter so that we reproduce the normalization of the empirical 
SK star formation law, which gives esF = 13Gyr~^. 

(5) Feedback: eefb, £mfb, Vsn, Esn, Psn- For the en- 
ergy driven wind model the mass loading factor is given by 
77 — 2eEFB SsN ??SN/T4ac; wMle for the momentum driven 
wind model r] = emfb Psn »?sn/Vosc. We adopt an energy 
per supernova of Esn = 10^^ erg, a momentum per super- 
nova of Psn = 3 X 10'* Mq kms~*, and a supernova rate of 
VSN = 8.3 X 10"'' per Solar mass of stars formed. For our en- 



ergy driven outflow model we adopt eefb = 0.25 (i.e. 25% of 
the SN energy ends up in the outflow) and for our momen- 
tum driven outflow model we adopt emfb = 1 (i.e. 100% of 
the SN momentum ends up in the outflow) . The no feedback 
model has eefb = emfb = 0. 



(6) Stellar populations and chemical enrichment: 
7?,, t/, Zhot, and the choice of initial mass function (IMF). 
We assume that the baryonic gas that is accreted by a dark 
matter halo has a Helium mass abundance of 0.25 and an ini- 
tial metallicity of Zhot = 0.002 ~ O.IZq. Contrary to DB09, 
where we adopted the instantaneous recycling approxima- 
tion (IRA), here we calculate the amount of recycled gas as 
a function of time by convolving the star formation history 
with the return fractions from the BC03 models. At low red- 
shifts this makes no difference to the final stellar masses, but 
at high redshifts the IRA over-estimates the return fraction 
and hence under-estimates the stellar masses. For simplic- 
ity, we still use the IRA for the chemical evolution, using a 
stellar yield of y = 0.02 appropriate for the Chabrier (2003) 
IMF adopted throughout. 



© 2009 RAS, MNRAS 000.[Tlf23l 



6 Button et al. 




9 10 11 9 10 11 9 10 11 



log,o Mstor [Mg] 

Figure 3. The SFR sequence for different feedback models: no feedback (red dots), momentum driven feedback (green dots) and energy 
driven (black dots). The slope, zero point and scatter of the SFR sequence is independent of feedback (the different colored dots trace 
the same relations). 



3 ORIGIN OF THE ZERO-POINT 

EVOLUTION OF THE SFR SEQUENCE 

3.1 Summary of observations 

The observational measurements of the SFR sequence from 
low to high redshifts that we use in this paper are sum- 
marized in Table [T] and Fig. [T] The slope of the SFR se- 
quence defined as dlog 5'Fi?/dlog Mstar is ~ 0.65 — 0.95. 
Shallower slopes tend to be measured at low redshifts (Elbaz 
et al. 2007; Noeske et al. 2007; Salim et al. 2007; Schimi- 
novich et al. 2007; Bothwell et al. 2009). Studies above 
z ~ 1 tend to find steeper slopes of around 0.9, in par- 
ticular the recent study by Pannella et al. (2009) finds a 
slope of 0.95 ± 0.07 at z ~ 2 for stellar masses in the range 
Mstar — lO^*^'^"^^" Af0. It is not clear if this implies an evo- 
lution of the slope, or is due to a systematic difference in 
SFR indicators or sample selection used at different red- 
shifts. The intrinsic scatter is somewhat uncertain, but is 
estimated to be < 0.3 dex (e.g. Noeske et al. 2007a; Elbaz 
et al. 2007; Daddi et al. 2007). The zero-point evolves by a 
factor of ~ 20 from redshift 2: ~ 2 to z ~ 0, but is roughly 
constant from z ~ 7 to z ~ 2. Note, though, that the high 
redshift (2 > 4) data have not been corrected for extinction 
as observations indicate that there is little or no extinction 
in these galaxies (Gonzalez et al. 2009). Nevertheless, we 
caution that the observed SFRs at 2 > 4 may have been 
underestimated. 



These relations assumed either a Salpeter IMF or a 
Chabrier (2003) IMF, as indicated in Table [U For con- 
sistency, we correct the stellar masses and SFRs to a 
Chabrier (2003) IMF by subtracting 0.25 dex (e.g. Gonza- 
lez et al. 2009). Given that the slope of the SFR sequence 
is close to unity, both Salpeter and Chabrier (2003) IMF's 
result in very similar zero-points to the SFR sequence. 



3.2 Comparison with observations 

Fig. [2] shows the SFR sequence at redshifts z = 0, 1, 2, 3, 4 
and 6 for our model with energy driven feedback (EFB). 
The red solid lines show power-law fits to the SFR sequence 
for stellar masses between 10^ ^ Mstar ^ 10^^ Mq, with 
slope (6), zero point (a) and scatter (cr) as indicated. The 
black error bars show the median SFR in stellar mass bins 
of width 0.2 dex. The solid black lines enclose 68% of the 
model data in each stellar mass bin, and indicate that the 
scatter in independent of mass. The model galaxies have 
been uniformly sampled in redshift zero halo mass (with 
10" ^ Mvir ^ IO^-^'^/i-^Mq). At high redshifts the tail 
of galaxies to high stellar masses is caused by scatter in 
the MAH. In particular, the highest mass galaxies at high 
redshifts (z > 2) are all in early forming (high concentration) 
haloes. 

The observed relations are shown by blue lines. Dashed 
lines show the mean, dotted lines show the la scatter, where 
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Figure 4. Evolution in the specific star formation rates for stellar masses of logjg Mstar/fAf©] = 9.3 (left panel), 10.0 (middle panel), 
and 10.7 (right panel). The observed evolution is indicated by symbols (see Fig.[T]for references). The horizontal error bars indicate the 
redshift range of the data points, the vertical error bars indicate the estimated intrinsic scatter in SSFR. Black lines correspond to our 
models: energy driven outflows (EFB, solid lines), momentum driven outflows (MFB, dashed lines), and no outflows (NFB, dotted lines). 
The evolution in the models is the same for all three feedback models, and closely follows the (1 + z)^-^^ scaling (red long-dashed lines) 
for the dark matter accretion rates at fixed halo mass expected in a ACDM cosmology. The models are in reasonable agreement with 
the data, although they seem to under predict the SSFRs at 2 ~ 2, and over predict those at 2 > 4. 



we have adopted a scatter of 0.3 dex. Note that the model 
scatter is significantly smaller than this (0.12 ± 0.01 dex), 
but is independent of mass and redshift. We do not consider 
this a serious problem for the model, given that the system- 
atic uncertainties in the data are likely to be significant (i.e., 
the observed scatter has to be considered an upper limit on 
the intrinsic scatter). In addition, it is likely that our model 
underestimates the true scatter due to our simplified treat- 
ment of the halo mass accretion history (i.e. the assumption 
that there is a one-to-one relation between MAH and halo 
concentration). For example cosmological hydrodynamical 
simulations find a scatter of about 0.3 dex in the gas ac- 
cretion rates in haloes of mass IO^^Mq at z — 2.5 (Dekel 
et al. 2009). Such a scatter in gas accretion rates would ex- 
plain the full observed scatter in the SFR sequence. 

The slope of the SFR sequence in our model is close to 
but less than unity, and evolves to slightly lower values at 
higher redshifts. In fact, the model slope seems to be some- 
what too steep compared to the data at z = 0. It remains to 
be seen whether this discrepancy holds up with better data 
coming available. 

Fig. [3] shows the SFR sequence for our three feed- 
back models: no feedback (red); momentum driven feedback 
(green); and energy driven feedback (black). This shows that 
the slope, zero point and scatter of the SFR sequence are all 
independent of the feedback model! This feature is remark- 
able because every other galaxy scaling relation depends 
quite significantly on the detailed treatment of feedback (e.g. 
van den Bosch 2002; DB09; Button 2009). The origin of this 
intriguing characteristic of the SFR sequence is discussed in 
^ below. 

Fig. [4] shows the evolution in the specific star formation 
rate (SFR/Mstar) at three stellar masses: Mstar = 2x10^ Mo 
(left-hand panel), Mstar = 1 x 10^° Mq (middle panel), and 
Mstar = 5 X 10^° M0 (right-hand panel). The data (described 
above) are shown as symbols, while the black lines indicate 



our models (solid line, EFB; dashed line, MFB; dotted line, 
NFB). Note that the models predict substantial evolution 
in the SSFRs, in rough overall agreement with the data, 
and that there is very little difference between the three 
different models (as also seen in Fig. [31). There is some in- 
dication that the models systematically under predict the 
SSFRs at 1 < z < 2, especially for the less massive galax- 
ies. This may signal a potential shortcoming of the mod- 
els, though it seems also possible that the observations have 
over-estimated the SSFRs due to sample selection effects or 
other systematics (see Dave 2008 for a discussion of system- 
atic effects in SFRs and stellar masses). 

The discrepancies between the SSFRs in our models 
and the data at 2: ~ 2 are of similar (but smaller) magni- 
tude and direction as found by Dave (2008), using cosmo- 
logical hydrodynamical simulations. Dave (2008) proposed 
that an evolving IMF could reconcile the data with theory. 
A key prediction of this evolving IMF model is that the dis- 
crepancy between the observed and predicted (with fixed 
IMF) SSFRs should increase monotonically with redshift. 
However, since the work by Dave (2008) a number of ob- 
servational studies have presented SSFRs at z ^ 4 (Daddi 
et al. 2009a; Stark et al. 2009; Gonzalez et al. 2009). In par- 
ticular, the study of Daddi et al. (2009a) measured the SFR 
sequence for B-band dropouts at 2; ~ 4, which were found 
to be only marginally higher than those measured at z ~ 2. 
Our models match this 2: ~ 4 data point remarkably well. 
This suggests that the discrepancy at 2: ~ 2 is not due to 
a variable IMF, at least not as in the model advocated by 
Dave (2008). 

Current observations indicate that the SSFRs are con- 
stant from redshifts 2 ~ 7 to 2: ~ 2. This is on stark con- 
trast with our models, which predict that the SSFRs de- 
crease monotonically with decreasing redshift. Even though 
the models reproduce the observed SSFRs at 2 ~ 4, at 2 = 7 
the discrepancy between the models and the observations is 
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Figure 5. Evolution in the specific baryon accretion rates for stellar masses of logjQ Afstar/[AfQ] = 9.3 (left panel), 10.0 (middle panel), 
and 10.7 (right panel). The observed evolution in specific star formation rates is given by the colored symbols (as in Fig.[T]l. The specific 
baryon accretion rates, Macc/A^star, into dark matter haloes are given by the black lines, and the specific baryon accretion rates onto 
the galaxies are given by red lines. The line types correspond to different feedback models: energy driven outfiows (EFB, solid lines), 
momentum driven outflows (MFB, dashed lines), and no outflows (NFB, dotted lines). The evolution in baryon accretion rates closely 
follows the evolution in specific star formation rates. However, the zero point depends on the feedback model, especially at lower stellar 
masses. 



a factor of ~ 4. As mentioned in ij3.1l these very high-redshift 
data have not been corrected for dust extinction, and thus 
if there is substantial extinction in very high redshift galax- 
ies (as there is in lower redshift galaxies), then this could 
reconcile observation with theory. 

3.3 What drives the evolution in the zero point of 
the SFR sequence? 

Having established that our models reproduce the decrease 
in SSFRs by a factor of ~ 20 from redshifts z ^ A to z = Q, 
we now address the question: what drives this evolution? 

3.3.1 Gas accretion rate? 

It has been pointed out by previous authors (e.g. Noeske 
et al. 2007b) that the decline of the SFRs at a given Matar 
is compatible with the cosmological decline in the accretion 
rate onto dark matter haloes, which scales roughly as ~ 
(l + z)'^'^^ at a fixed halo mass (Birnboim, Dekel, & Neistein 
2007). Note that the actual slope of the specific accretion 
rate varies from 2 at z = to 2.5 at high redshift (Neistein 
& Dekel 2008). 

Fig. [5] shows the evolution of the baryon accretion 
rate onto the halo (Mace) vs. redshift (black lines) and the 
evolution of the cold gas accretion rate onto the galaxy 
(Mcooi) (red lines). As expected from Figs. IBl] fc lB21 these 
two curves only differ significantly at high stellar masses 
(Mstar > 5 X 10^" M0), where the cooling time becomes 
larger than the free fall time. We note that cold streams 
are likely to eliminate this difference. For comparison, the 
observed evolution in star formation rates are shown as col- 
ored symbols (as in Fig. |4]). The main outcome is that the 
evolution in the baryon/cold gas accretion rate closely fol- 
lows the evolution in the star formation rate (both for the 
models as well as in the data). This confirms that a plausible 



explanation for the decline in the SSFRs of galaxies is sim- 
ply a decline in the specific gas accretion rate. However, the 
zero point is dependent on the feedback model, especially 
for lower stellar masses. This is because feedback is more 
efhcient at removing baryons from lower mass haloes. For 
a given stellar mass, the halo mass will be higher for mod- 
els with stronger feedback. Since higher mass haloes have 
higher baryon accretion rates (in absolute terms), the mod- 
els with stronger feedback have higher specific baryon accre- 
tion rates. 

The upper panels of Fig. [S]show the SFR vs. the cold 
gas accretion rate onto the galaxy, A/cooi, for our three feed- 
back models. For the NFB model without outflows, there is 
an almost 1:1 correspondence between the SFR and inflow 
rate. However, for models with outflows (MFB and EFB), 
the relation between SFR and inflow rate has a slope greater 
than unity, with a steeper slope for the EFB model than for 
the MFB model. The differences between the three models 
are largest at low accretion rates, which correspond to low 
mass haloes, and hence where the speciflc outflow rates are 
highest. This suggests that the net inflow rate might corre- 
late better with the SFR. This is indeed the case, as shown 
in the lower panels of Fig. [B] Thus we conclude that galaxy 
disks are in a steady state in which the star formation rate 
is roughly equal to the net gas inflow rate. This conclusion 
can also be derived by solving the differential equation for 
change in gas mass with time (Bouche et al. 2009). 



3.3.2 Gas masses or gas densities? 

We have shown that the decline in net gas accretion rate 
can explain why SSFRs decline with time, but this does not 
explain how the decline in gas accretion rate is converted 
into a decline in the SSFRs. Recall that, as discussed in the 
introduction, the evolution of the star formation rate at a 
given stellar mass is expected to depend on the change in 
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Figure 6. The SFR as function of the rate at which cold gas is accreted by the galaxy (Mcooi; upper panels) and as a function of the net 
inflow rate of cold gas onto the galaxy, given by Mcool — Mout (lower panels). Results are shown for model galaxies with stellar masses 
Matar > 10^ Mq and for our three feedback models: energy driven outflows (EFB, black dots), momentum driven outflows (MFB, green 
dots), and no outflows (NFB, red dots). For the NFB model the SFR is proportional (but not identically equal to) the cold gas inflow 
rate. For models with outflows, the SFR is not directly proportional to the gas inflow rate. For these models the SFR is better correlated 
with the net gas inflow rate. This indicates that our model galaxies have settled in a steady state between the rate of inflowing gas and 
the rate of star formation. 



gas mass and/or gas density at that stellar mass. A factor 
20 decrease in SFR needs a factor 20 decrease in gas mass, 
or a factor 400 decrease in gas surface density (assuming a 
SK star formation law with n — 1.5) or factor of 20 decrease 
in gas surface density (for a SK law with n = 2). Note that 
while the global SK law has a slope of n ~ 1.5, for spiral 
galaxies the spatially resolved SK law typically has n ~ 2 
(Bigiel et al. 2008, and references therein) as thus it may be 
the more appropriate value to use to estimate the required 
evolution in gas disk density. 

Fig. [7] shows the evolution in the cold gas to stellar 
mass ratio, Mgas/Mstar, where cold gas is both atomic and 
molecular. In all models the evolution is weak, especially 
for the two models with outflows (solid and dashed lines). 
The outflow models also reproduce the cold gas fractions 
at redshift z = 0, whereas the no feedback model (dotted 
lines) systematically under-predicts the cold gas masses at 
all stellar masses. At high redshifts the observations are not 
direct and thus carry significant uncertainties. At z ~ 2, Erb 
et al. (2006) derive gas masses from the SK star formation 



law. In addition to the systematic uncertainties (such as un- 
certainties in gas disk sizes, star formation rates, and stellar 
masses) inherent in these measurements there are concerns 
about sample selection. Since the Erb et al. (2006) sample 
is UV selected, it may be biased towards gas rich, high star 
formation rate galaxies, especially at lower stellar masses. 
Given these uncertainties, we consider our models are not 
necessarily in disagreement with the observations. 

At 2 ~ 1.5, Daddi et al. (2009b) detected the C0[2- 
1] line for 6 star forming galaxies, and inferred a mean gas 
mass fraction of 57±6%. These gas fractions are higher than 
predicted by our model. However, there is a large uncertainty 
(a factor of ~ 4) in the conversion factor, aco , between CO 
luminosity and molecular gas mass. Daddi et al. (2009b) 
attempt to infer aco using dynamical models, but these 
involve many assumptions, including the (unknown) dark 
matter fraction, and thus do not robustly constrain qco- 
We note that even with the high qco advocated by Daddi 
et al. (2009b) this only results in a factor of ~ 3 decrease 
in gas masses between z ~ 1.5 and z = 0. A factor of ~ 10 
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Figure 7. Evolution in tlie ratio between cold gas and stellar mass for stellar masses of logjQ Mstar/[MQ] = 9.3 (left panel), 10.0 
(middle panel), and 10.7 (right panel). The models are given by black lines as indicated. The observed gas-to-star ratios are given by 
blue symbols as indicated. At low redshift the data reflects direct measurements of atomic and molecular gas (from DB09 using data 
from Garnett 2002), and are in good agreement with the energy and momentum driven wind models. At z = 2 the data reflect gas 
mass fraction derived indirectly using the Schimdt-Kennicutt star formation law (Erb et al. 2006, E06). At 2 = 1.5 the gas fractions are 
based on CO measurements and dynamical models (Daddi et al. 2009b, D09). The small (blue) error bar gives the nominal uncertainty 
while the large (red) error bar shows the uncertainty in the conversion factor from CO luminosity to molecular gas mass. At z = 1.2 and 
2 = 2.3 the gas fractions are based on CO measurements (Tacconi et al. 2010, TIO). The evolution in the model gas fractions is weak, 
especially for the two models with outflows. We conclude that the evolution in the SFR sequence is not dominated by evolution in the 
cold gas mass fractions. 
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Figure 8. Evolution in the effective surface density of the cold gas (defined as the average density of cold gas within the cold gas half 
mass radius) for stellar masses of logj^Q Mstar/[MQ] = 9.3 (left panel), 10.0 (middle panel), and 10.7 (right panel). The models are given 
by black lines, as indicated. There is strong evolution in the effective gas densities for all three models and for all three stellar masses. 
The red dashed lines show the evolution of the halo surface density at fixed halo mass. Given the weak evolution in the gas- to-stellar 
mass ratios (Fig.[7]l, it is this evolution in the gas surface densities that provides the mechanism for the evolution in gas accretion rates 
to turn into an evolution in the SFR sequence. 



decrease in gas masses is required for change in gas fractions 
to explain the factor of ~ 10 decrease in SSFR between z = Q 
and z ~ 1.5. 

More recent observations (Tacconi et aL 2010) have de- 
tected the CO [3-2] hne for 19 massive (Mstar ~ 5 x 10^° - 
3 X lO^^M©) star forming galaxies at redshifts z ~ 1.2 and 
z ~ 2.3, finding mean molecular gas fractions of 0.34 and 
0.44, respectively. These measurements fall between those 
of Erb et al. (2006) and Daddi et al. (2009b), and are thus 
also consistent with our models. 



We conclude that the current observations of gas frac- 
tions are consistent with weak or no evolution between z ~ 2 
and z = Q,as our models predict. But obviously this needs to 
be confirmed with direct observations of atomic and molec- 
ular gas masses at high redshifts, for larger samples of star 
forming galaxies and over a wider range of stellar masses. 

Fig. [8] shows the evolution in the effective gas surface 
density, which we define as the mean gas density within 
the radius enclosing half of the cold gas mass: Eso.gas = 
Mgas/(27ri?|o gas), where J?5o,gas is the half mass radius of 
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Figure 9. Evolution in the cold gas-to-stellar mass ratios for the EFB model, where the line types correspond to the atomic (dotted, 
light grey), molecular (dashed, red) and total (atomic plus molecular; solid, dark grey). The three panels correspond to stellar masses 
of logxo Mstar/[A^Q] = 9.3 (left panel), 10.0 (middle panels), and 10.7 (right panels). At fixed stellar mass, the molecular gas masses 
increase by an order of magnitude (independent of the details of gas accretion and feedback) , while the atomic gas masses decrease by 
a similar amount. Since molecular gas is the fuel for star formation, and the SFR ~ M^^l^ the evolution in the molecular gas masses at 
fixed stellar mass translates directly into evolution in the specific star formation rates. 



the gas. In our models the gas density shows very strong 
evolution: a factor of ~ 10 from z — 2 to z — 0, and a 
factor of ~ 100 from z = 4 to z = 0. It is thus the increase 
in the density of the gas that provides the mechanism for 
turning the higher net gas accretion rates at high redshifts 
into higher star formation rates, at least in our models. 



3.3.3 Why are high redshift disks denser? 

Obviously this poses the question why gas disks at higher 
redshifts have higher surface densities. A crude approxima- 
tion to the evolution of disk sizes, Rd, can be made by assum- 
ing that 1) galaxies have fiat rotation curves and 2) that the 
disk spin parameter, Ad, is independent of redshift. In that 
case, the disk size scales with the size of the dark matter 
halo according to Rd{z) ^ Ad-Rvir(z) (Mo, Mao, & White 
1998). The size of the halo scales with mass and redshift 
as R^i,{z) ~ 7i"(2)-2/3Avir(2)"^/^Mvir(z)'''^ Here H{z) is 
the Hubble parameter and Avir(z) is the overdensity of the 
halo with respect to the critical density of the Universe. 
Thus if the size and mass of the gas disk are proportional 
to those of the halo, then the density of the gas disk scales 
as T,{z) ~ H{z)'^/^A^ir{zf^^ at fixed halo mass. This scal- 
ing is shown (with arbitrary normalization) as the red long- 
dashed lines in Fig.[S] Our model roughly follows this scaling 
at 2 < 1 — 2.5 (depending on stellar mass). For higher red- 
shifts, however, our model predicts an evolution in density 
that is significantly stronger. This indicates that there are 
other effects in our model that drive the evolution of gas den- 
sities besides simply halo density evolution. As discussed in 
Appendix [C] the predominant additional effect is a decrease 
in the disk specific angular momentum at high redshifts. 

It should be noted that in our models the gas densities 
are higher at higher redshift because the half mass radii are 
smaller. An equivalent effect could in principle be achieved 
if gas disks at high redshift were considerably more clumpy 
than low redshift gas disks of the same stellar mass. The 



dumpiness could result in a larger fraction of the gas disk 
being above the threshold for efficient molecular gas forma- 
tion, and thus result in higher molecular gas masses and 
higher star formation rates, even though the sizes of the gas 
disks may not be significantly smaller than at lower red- 
shifts. This may be relevant given recent observational and 
theoretical evidence for clumpy disks at high redshifts (e.g. 
Elmegreen et al. 2007; 2009; Agertz, Teyssier, & Moore 2009; 
Ceverino, Dekel, & Bournaud 2009; Dekel, Sari, & Ceverino 
2009). 



3.3.4 Molecular vs. atomic gas 

We have shown that while the gas fractions show only weak 
evolution with redshift, gas densities are much higher at 
higher redshift. Hence, one would expect that the molecular 
gas fractions will also be higher at higher redshift. Fig. [5] 
shows the ratios of total, atomic and molecular gas to stel- 
lar mass for the EFB model as functions of redshift. As 
expected, at a fixed stellar mass, the molecular gas masses 
increase by an order of magnitude from low to high redshift, 
whereas the atomic gas masses decrease by a corresponding 
amount. Recent observations of molecular gas masses in star 
forming galaxies at redshifts z ~ 1.2 — 2.3 support this result 
(Daddi et al. (2009b); Tacconi et al. 2010). Given that the 
Schmidt star formation law for molecular gas has a slope 
~ 1 (e.g. Bigiel et al. 2008), this implies that the SFRs are 
roughly proportional to molecular gas mass (Eq.4). Thus the 
increase in molecular gas fraction provides the fuel for the 
higher star formation rates at higher redshifts. 

Fig. [TO] shows the evolution in the effective gas surface 
density for the EFB model. The solid line corresponds to 
the total gas as in Fig. [S] The dashed lines show the molec- 
ular gas density, while the dotted lines show the atomic gas 
density. The molecular gas density shows a strong evolu- 
tion with redshift, whereas the atomic gas density remains 
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Figure 10. Evolution in the effective gas surface density for the EFB model, where the line types correspond to the atomic (dotted), 
molecular (dashed), and total (atomic plus molecular; solid) gas. The three panels correspond to three different stellar masses, as 
indicated. Note that the surface densities of atomic gas are virtually constant (at ^ IOMqpc"^) while the surface densities of molecular 
gas evolve strongly with redshift. This reflects the fact that at surface densities in excess of ~ lOM0pc~^, the midplane pressure in the 
disk is sufficiently high to assure an efficient conversion of atomic to molecular gas (see Appendix A). 



roughly constant at Eso.gas — IOMqpc^ all the way from 
2; ~ 6 to 2 = 0. 

In the near future the Atacama Large Millime- 
ter /submillimeter Array (ALMA) will have the sensitivity to 
test our theoretical prediction that the molecular gas den- 
sities and molecular gas fractions of disk galaxies increase 
strongly with redshift. A test of our predictions for the evo- 
lution of the mass fractions and densities of atomic gas, how- 
ever, most likely will have to wait for the Square Kilometer 
Array (SKA), which is currently scheduled to start operation 
no earlier than 2018. In the meantime, measurements of HI 
column densities from damped Lya systems (DLAs) provide 
a means of measuring the atomic gas content of high redshift 
galaxies. The analysis of 738 DLAs by Prochaska & Wolfe 
(2009) suggests that HI disks have not evolved significantly 
since 2; ~ 2. Although in agreement with our model pre- 
dictions, a more detailed comparison between the observed 
DLA column densities and our model predictions is needed 
to assess whether they are indeed consistent with each other. 

In summary, the increase in net gas accretion rate from 
low to high redshift provides the means for the SFRs to in- 
crease at a fixed stellar mass. The mechanism for converting 
this elevated gas accretion rate into an elevated star forma- 
tion rate is an increase in gas density, not an increase in 
gas mass. Increased gas density, mainly due to the fact that 
haloes are denser at higher redshifts, results in increased 
molecular gas fractions, which provides the increased fuel 
supply for the higher SFRs at high redshifts. 



4 ORIGIN OF SCATTER IN THE SFR 
SEQUENCE 

In this section we address the question of why the scatter 
in the SFR sequence is so small. In our model there are 3 
sources of scatter: halo concentration, c, halo spin. A, and 
halo AMD, a. We assume the scatter in these three pa- 
rameters are independent of each other, and independent of 
redshift. 



Fig. [TT] shows the SFR sequence (expressed in terms of 
SSFR vs. Matar) at redsMfts z = 0, 1, 2, and 4 for the EFB 
model. Each model has only one source of scatter: halo con- 
centration (0.11 dex) halo spin (0.22 dex), halo AMD (0.11 
dex). This shows that the scatter in the SFR sequence, at all 
redshifts is dominated by the scatter in halo concentration, 
and rather remarkably, shows only a very weak dependence 
on halo spin or AMD. 



4.1 Why does the scatter in the SFR sequence 
depend on halo concentration? 

In our model the concentration of the dark halo is directly 
coupled to the mass accretion history (MAH) of the dark 
matter halo (see DB09 for details). In order to assess the 
relative importance of scatter in MAH vs. scatter in halo 
concentration, we have also run models in which the con- 
centration is decoupled from the MAH. These show that 
scatter in the halo concentration parameter contributes only 
a small amount to the scatter in the SFR sequence. Thus it 
is mainly the variations in the MAHs that dominates the 
scatter in the SFR sequence, at least in our models. Higher 
concentration halos have MAHs shifted to earlier times (i.e., 
these haloes assemble their mass earlier), and vice versa for 
lower concentration halos. An earlier MAH results in higher 
baryon accretion rates at high redshifts, which in turn im- 
plies higher SFRs, but also a larger stellar mass. 

In our model the scatter in MAH is determined by scat- 
ter in halo concentration, which is 0.14 dex of all haloes 
(Bullock et al. 2001a) and just 0.11 dex for relaxed haloes 
(Wechsler et al. 2002; Maccio et al. 2007). Since we adopt 
a smooth MAH, and we assume a one-to-one correspon- 
dence between halo concentration and formation redshift, 
it is likely that our model underestimates the true scatter 
in baryon accretion rates onto central galaxies. Determining 
the true scatter requires cosmological hydrodynamical sim- 
ulations. The simulations used in Dekel et al. (2009) find 
a scatter of about 0.3 dex in the gas accretion rates onto 
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Figure 11. The impact of different sources of scatter on the scatter in the SFR sequence for the EFB modeh halo concentration (left 
panel); halo spin parameter (middle panel); halo AMD parameter (right panel). The model galaxies at redshifts z = 0,1,2, & 4 are 
indicated by grey points, the black error bars show the error on the median SSFR in Mstar bins of 0.2 dex width, and the black lines 
enclose 68.3% of the model galaxies. The average scatter in SSFR over the four redshifts is given in each panel. The scatter in SSFRs 
is clearly dominated by scatter in the halo concentrations, while scatter in the spin parameter or AMD parameter result in negligible 
scatter in SSFRs. 



galaxies in haloes of mass 10^^ Mq at z — 2.5. Assuming 
this scatter translates linearly into scatter in SSFR (as in 
our models) this would reproduce the full observed scatter 
in the SFR sequence. Further studies are needed to quantify 
the mass and redshift dependence of this scatter. 

The upper panels of Fig. [T^] show the effect of varying 
halo concentration c (and thus also the MAH), halo spin 
parameter, A, and AMD, a, on the SFR sequence. The solid 
black circle shows a model with median c, A and a, and 
with energy driven feedback with ere = 0.25. The open 
blue circles show models with ±2a variation in c. The di- 
rection the model moves for increasing c is given by the 
blue arrows. This shows that galaxies that form in haloes 
with higher concentrations (i.e. earlier MAHs) lie below the 
SFR sequence (i.e. lower SFR at fixed Mstar), while galaxies 
that form in haloes with lower concentrations lie above the 
SFR sequence (i.e. higher SFR at fixed Mstar). This implies 
that the normalization of the SFR sequence depends on the 
cosmological parameters, in particular ag, with higher erg 
resulting in lower SFR at a given i\fstar (see § I5.3|l 



4.2 Why is the scatter in the SFR sequence 
independent of halo spin? 

The red open squares in Fig.[T2]show models with ±2a vari- 
ation in A, while the green open triangles show models with 
±2cr variation in a. As for the concentrations, the direction 
the model moves for increasing A and a is indicated by the 
red and green arrows, respectively. This shows that both 
SFR and Mstar depend on the halo angular momentum pa- 
rameters, and spin in particular, but that the changes move 
galaxies along the SFR sequence: lower A and a both result 
in higher SFR and higher Mstar . 

The middle and lower panels of Fig. [12] show the effect 
of variation in c. A, and a, on the relation between stel- 
lar mass and the gas-to-stellar mass ratio (middle panels) 
and between stellar mass and gas density (lower panels). 



The scatter in these two relations is dominated by A, with 
c a sub-dominate component, and a a marginal component. 
Changes in a (AMD) only have a very small effect, which 
explains why in many cases you can only see one green tri- 
angle (that for low values of a). Smaller values for A result 
in higher gas densities, which result in more efficient star 
formation, and hence higher SFRs, higher Mstar, and lower 
Mgas . This explains why the SFR sequence is invariant under 
changes in A. 



4.3 Independence of the SFR sequence to disk 
size 

In our model the scatter in disk sizes, at fixed stellar mass, is 
dominated by scatter in the spin parameter. Scatter in halo 
concentration (which is coupled to scatter in MAH) also 
contributes, but to a much lesser extent (see for example 
Dutton et al. 2007, Fig. 7). Thus we would expect that the 
SFR sequence is also independent of disk size. In Fig. [TS] 
we show that this is indeed the case. This figure shows the 
scatter in SSFR at fixed stellar mass vs. the scatter in disk 
size at fixed stellar mass, for galaxies in the EFB model 
with stellar masses between 9.0 < logjQ Mstar < 10.5. We 
get similar results, i.e. no correlation, for our other feedback 
models. In the upper panels the disk size is the stellar half 
mass radius, while in the lower panels the disk size is the 
cold gas half mass radius. The left panels show results at 
redshift z = 0, the middle panels z — 1 and the right panels 
z = 2. 

We now return to the line of reasoning from the intro- 
duction, that suggested that scatter in disk size would be 
expected to result in scatter in the SFR sequence, and ex- 
plain why it is false (at least in our model). Recall that for 
a SK star formation law with n = 1.5 the SFR depends lin- 
early on gas mass and gas disk size. Thus for fixed gas mass, 
one expects that galaxies with smaller gas disks should have 
higher SFRs. If stellar mass is proportional to gas mass, and 
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Figure 13. Correlation between the scatter in specific star formation rate (SSFR) at fixed stellar mass to scatter in disk size at fixed 
stellar mass in the EFB model at redshifts z = 0, 1, & 2. The disk sizes arc stellar half mass radii (upper panels) and cold gas half mass 
radii (lower panels). This shows that the scatter in the SFR sequence is independent of disk size (both stellar and gaseous), as expected 
from the independence of the SFR sequence to spin parameter. 



stellar disk sizes are proportional to gas disk sizes, then one 
would expect that at fixed stellar mass, scatter in stellar disk 
sizes to result in scatter in SFRs. The reason that this is not 
the case in our models is that the scatter in the relation 
between gas mass and stellar mass depends on disk size. At 
a fixed stellar mass, galaxies with smaller disks have lower 
gas fractions, and galaxies with larger disks have higher gas 
fractions. These two effects cancel out, leaving the SFR un- 
changed. 

Another way to think about this is that at fixed halo 
mass there is a spread in disk sizes. Galaxies with smaller 
disks do have higher SFR, but at the same time they also 
have higher stellar masses (and lower gas fractions). Thus 
scatter in disk sizes creates a SFR sequence at fixed halo 
mass. 

The lack of a correlation we find between scatter in 
SSFR and scatter in disk size (both stellar and gas) at fixed 
stellar mass should be testable using existing and future 
data. Indeed, correlations between SSFR and stellar sur- 
face density have been studied by Franx et al. (2008) and 
Williams et al. (2009). These authors have shown that there 
is a strong correlation between SSFR and stellar surface den- 
sity, such that galaxies with higher stellar densities have 



lower SSFR. Williams et al. (2009) also studied the galaxy 
half fight radius - stellar mass relation as a function of SSFR 
from redshifts z = to z = 2. At a fixed stellar mass, there 
is a clear trend for smaller galaxies to have lower SSFR. 
These observations may seem at odds with the predictions 
from our models. However, our predictions are for galax- 
ies on the SFR sequence, and not for galaxies in general. 
Thus it is plausible that the stellar surface density is impor- 
tant for determining whether a galaxy is on or off the SFR 
sequence. But that for galaxies on the SFR sequence, the 
stellar density does not correlate with the SSFR, as in our 
models. Another caveat is that the sizes measured by Franx 
et al. (2008) and Wilhams et al. (2009) are global half light 
sizes, whereas our model sizes are disk sizes. Global sizes do 
not distinguish between bulges and disks. Since disks tend 
to be larger than bulges as well as having higher SSFRs, a 
variation of bulge fraction at fixed stellar mass will naturally 
result in a correlation between galaxy half hght radius and 
SSFR, even if no such correlation exists between disk size 
and SSFR. Thus averaging bulges and disks together may 
be hiding some important clues to the strength and nature 
of the correlation between SSFR and stellar surface density. 



© 2009 RAS, MNRAS 000.[Tlf23l 



16 Button et al. 



a: 



o 



- 


1 ' 


1 ' 1 ' . 
e„= 1 3 _ 

or 






H-2 - - - 
x2 - 




z = 2 






z=1 




- z=0 - 




1 . 


1,1," 



1 - 



- = 



T 1 1 1 1 1 ■- 

n.=0.022 — _ 
" -2-- 
x2 



1 - 



- 



T — ' — r 

(78 = 0.8- 

-0.2- 
,-1-0.2 ■ 



9 10 11 
'oQio M,.or [Mo] 



9 10 11 
'09,0 M,.o, [M©] 



9 10 11 
"09.0 [M^] 



Figure 14. Dependence of the specific star formation rate (SSFR) - stellar mass relation on the star formation efficiency parameter, 
egp (left panel); cosmological baryon density, Q^, (middle panel); and the cosmological power spectrum amplitude, erg (right panel). The 
relation due to our adopted EFB model is given by the solid black lines. Relations due to lower /higher values of these parameters are 
given by dashed green and dotted red lines, respectively. The SSFR - Afstar relation is independent to changes in the star formation 
efficiency parameter (at the factor of 2 level), and also largely independent of the cosmological baryon density. The zero point of the 
SSFR- A/star relation is dependent on erg, as this parameter controls the timescale of the halo mass accretion history. 



Future observational studies will be able to address these 
issues. 



5 ROBUSTNESS OF THE SFR SEQUENCE 

5.1 Dependence on feedback 

The black open pentagons in the upper panels of Fig. [12] 
show the effect of varying the energy feedback efficiency, 
eefb, from to 1 on the SFR sequence (eEPB increases 
from right to left as indicated by the black arrow). Varying 
ebfb moves galaxies along the SFR sequence, in the sense 
that models with stronger feedback have both lower stellar 
masses and lower SFRs, but the same SSFR. Increased feed- 
back efficiency results in lower gas densities (lower panels), 
and hence less efficient star formation (upper panels), and 
higher gas-to-stellar mass ratios (middle panels). 

The fact that feedback results in higher cold gas frac- 
tions might seem counter-intuitive, as feedback blows cold 
gas out of the system which lowers the cold gas to halo mass 
fraction. However, the outflows in our model occur over an 
extended period of time, and thus the galaxy has time to 
adjust to them. Outflows reduce the baryonic masses, while 
at the same time increasing the specific angular momentum 
of the remaining baryons (Mailer & Dekel 2002; Button & 
van den Bosch 2009). These two effects cause the density of 
the gas disk to decrease, resulting in less efficient star forma- 
tion and hence higher cold gas fractions for the remaining 
baryons. In fact, it is interesting to note that varying the 
feedback efficiency moves galaxies roughly along a line of 
constant gas mass (red long-dashed line in Fig. I12[). 

5.2 Dependence on star formation efficiency 

One of the uncertainties in our star formation law is the 
normalization, which depends on amongst other things, the 
stellar IMF. The right panel of Fig. [H shows that the SFR 



sequence is independent to changes in the star formation 
efficiency parameter (esp) by a factor of 2. Furthermore, for 
individual galaxies, SFR and A/star are insensitive to esp. 
This is because changes (e.g. an increase) in esp are offset 
by changes (e.g. a decrease) in gas masses and gas densities, 
i.e. star formation is self-regulating. Thus differences in IMF, 
such as between Salpeter and Chabrier (2002), do not impact 
on the SFR sequence in our models. 



5.3 Dependence on cosmological parameters 

As discussed in tj4. li the SFR sequence depends on the MAH 
of the baryons, which is cosmology dependent. The MAH 
can be described to lowest order by a normalization, and a 
formation redshift. 

We simulate variation in the normalization by changing 
the cosmological baryon fraction by plus or minus a factor 
of 2. As shown in the middle panel of Fig.[T4]this has no ef- 
fect on the normalization of the SSFRs. The only effect is to 
change the mass scale at which cooling becomes inefficient, 
which occurs because changing the baryon fraction changes 
the density of the hot halo gas, and hence the cooling effi- 
ciency. 

We simulate variation in the formation redshift of the 
haloes by varying the normalization of the power spectrum, 
(Tg. Higher as results in earlier forming haloes, and hence an 
earlier MAH. The right panel of Fig. ll4l shows the impact on 
the SFR sequence of varying erg from 0.6 to 1.0. The SSFRs 
are higher for lower as, and lower for higher as- 

Thus the zero point and evolution of the SFR sequence 
appears to depend only on the time dependence of the gas 
mass accretion history. This is a result of galaxies being in a 
quasi steady state between gas accretion rate and star for- 
mation rate, which in turn is the result of the self-regulating 
nature of star formation in galaxy disks. This suggests that 
in order to change the slope or time evolution of the SFR 
sequence relation, it is necessary to change the gas accretion 
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history, rather than the details of how gas is converted into 
stars. 



5.4 Is the SFR sequence special? 

We have shown that the SFR sequence in our models is inde- 
pendent of the feedback efficiency, and has small scatter. We 
now consider relations between the SFR and other masses, 
to see if the SFR sequence is special. 

Fig. [TSl shows the SFR - M relations at redshift 2 = 0, 
where M is the virial mass (Mvir), galaxy mass (Afgai — 

Mstar + Mgas), StcUar mass (Mstar), cold gas mass (Afgas = 

Matom +Mnioi), cold atomic gas (Matom), and cold molecu- 
lar gas (Mmoi)- In each panel results are shown for all three 
models; NFB, EFB and MFB (as indicated). As we have al- 
ready seen, the SFR sequence (upper right panel) has small 
scatter (0.11 dex in Mstar at SFR = 0.1 to 10 yr"^) and 
a slope that is independent of mass and feedback. The SFR 
- Mgai relation (upper middle panel) has a slightly larger 
scatter (0.12 dex), which increases at the low mass end, and 
has a slope that is weakly dependent on mass and feedback. 
The SFR - Mvir relation (upper left panel) has a slightly 
smaller scatter (0.10 dex) than the SFR sequence, but the 
slope is strongly dependent on mass and the feedback model. 



The SFR - Mgas relation (lower left panel) has a large scat- 
ter (0.26 dex), and a slope that depends on the feedback 
model. The SFR - Matom relation (lower middle panel) has 
the largest scatter (0.32 dex) and also has a slope that de- 
pends on feedback. Finally, the SFR - Mmoi relation (lower 
right panel) has the smallest scatter (0.06 dex) and is inde- 
pendent of feedback. 

The origin of the tight relation between SFR and Afmoi 
is due to our effective star formation law, Esfr oc Smoi, 
with a Schmidt-law index of n ~ 1. Thus from Eq.(4), the 
integrated SFR of a galaxy (which reflects its stellar mass) 
is proportional to the total molecular gas mass. 



6 SUMMARY 

We have used a disk galaxy evolution model to investigate 
the origin of the zero point evolution and small scatter of the 
star formation rate sequence — the tight correlation between 
galaxy star formation rate (SFR) and stellar mass (Mstar). 
In our model the mass accretion histories of the dark matter 
haloes are smooth (i.e. there are no mergers), gas accretes 
onto the central galaxy via "hot mode" accretion (i.e. gas 
that enters the halo shock heats to the virial temperature, 
cools radiatively, and then accretes onto the central galaxy 
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in a free fall time), the cold gas density profile is determined 
by detailed conservation of specific angular momentum, the 
radially dependent molecular gas fraction is taken to depend 
on the mid-plane pressure in the disk, and star formation 
occurs at a rate that is governed by the local surface density 
of the molecular gas. We summarize our results as follows: 

• The slope, zero point and scatter of the SFR sequence 
are independent of our feedback model. This is because 
feedback effects the SFRs and Mstar such that galaxies are 
shifted along the relation, rather than perpendicular to it. 

• The zero point of the SFR sequence in our models 
evolves strongly, and in good agreement with observations 
between z ~ 4 to a = 0. However, the models under pre- 
dict the SFRs at 2; ~ 2, in qualitative agreement with the 
hydrodynamical simulations of Dave (2008) , who argued for 
evolution in the IMF in order to reconcile his simulations 
with the data. However, our model fits the data at a ~ 4, 
disfavoring such an evolving IMF model. 

• In our model the evolution in the zero point of the SFR 
sequence closely follows the evolution in the gas accretion 
rate. However, at a fixed stellar mass the absolute value of 
specific gas accretion rate depends strongly on the feedback 
model. There is a better correlation between the star for- 
mation rate and the net gas inflow rate (inflow - outflow). 
The high (central) densities of ACDM haloes and the rela- 
tively low values of the spin parameter aspire to assure that 
the majority of (centrifugally supported) disks always attain 
sufficiently high gas densities for star formation to be effi- 
cient, and thus for a steady state to exist in which the star 
formation rate follows the net inflow rate of cold gas. 

• Our models predict that the SSFR should decrease 
monotonically from redshifts 2; ~ 7 to z ~ 2, whereas the 
current observations favor no evolution over this redshift 
range. At ^ ~ 7 the discrepancy between models and ob- 
servation is a factor of ~ 4. Thus, if the observations axe 
correct, this indicates that the star formation rate no longer 
follows the gas accretion rate in high redshift {z > 2) galax- 
ies, as we infer that it does in low redshift {z < 2) galaxies. 

• At a fixed stellar mass, gas masses evolve only very 
weakly in our model, whereas gas densities evolve strongly. 
Thus the higher SSFRs in higher redshift galaxies is due to 
higher gas densities, and perhaps surprisingly, not due to 
higher total cold gas masses (as is often claimed e.g. Daddi 
et al. 2009b). The higher gas densities are driven by the 
higher densities of ACDM haloes at higher redshifts coupled 
to a decrease in disk specific angular momentum at higher 
redshifts. 

• The increased gas densities at higher redshifts result in 
higher molecular gas masses. At a fixed stellar mass, the 

molecular gas masses are a factor of ~ 10 higher at z ~ 3 
compared to 2 = 0. This thus provides the increased fuel for 
the higher SSFRs in high redshift galaxies. 

• At fixed stellar mass, the molecular gas densities in our 
model increase strongly with increasing redshift, whereas 
the atomic gas densities are roughly independent of redshift, 
with an efl^ective HI surface density of Sso.hi — 10 Mq pc~^. 
This is basically a reflection of the fact that at HI surface 
densities higher than this, the atomic gas is efficiently con- 
verted into molecular gas (see also Martin & Kennicutt 2001; 
Wong & Blitz 2002; Bigiel et al. 2008). 

• The scatter in SFR at fixed Mstar in our models is 



~ 0.12 ± 0.01 dex. The source of this scatter is the vari- 
ance in mass accretion histories for haloes of a given virial 
mass, which in our model is coupled to the halo concen- 
tration. Scatter in halo spin contributes negligibly, because 
spin scatters galax:ies along the SFR sequence. Note that 
the 0.12 ± 0.01 dex scatter in our models is less than the 
observed scatter 0.3 dex), as it should be given that the 
observed scatter is likely to be dominated by observational 
uncertainties. However, it is likely that our model underesti- 
mates the scatter due to our simplified treatment of the halo 
mass accretion history. It remains to be determined whether 
the true intrinsic scatter is indeed as small as predicted by 
our models, or whether additional sources of intrinsic scatter 
are required. 

• Relations between SFR and other masses (halo virial 
mass, galaxy mass, gas mass, atomic gas mass) exist, but 
have either larger scatter, or a slope that is dependent on 
mass and/or the feedback model. The tightest correlation is 
between SFR and molecular mass, which is a manifestation 
of the approximately linear relation between SFR surface 
density and molecular gas surface density which is built into 
our models. 

While our model is conceptually very simple, it re- 
produces the main features of the observed SFR sequence: 
strong zero point evolution, small scatter, and a slope that 
is close to unity. In detail, though, there seem to be three 
potential failures: (i) the model slope appears somewhat too 
steep at « = 0, (ii) the model underpredicts the SSFRs at 
« ~ 2, and (iii) the model overpredicts the SSFRs at z > 4. 
Better data is required in order to asses the severity of these 
potential shortcomings. 

A key prediction of our model that will soon be testable 
observationally (with ALMA) is that, at fixed stellar mass, 
the molecular gas densities and molecular gas masses should 
be substantially (by at least an order of magnitude) higher 
at z > 2 than they are for present-day disk galaxies. Results 
for small numbers of disk galEixies at 2 ~ 1 — 2 lend support 
to our predictions (Daddi et al. 2009b; Tacconi et al. 2010), 
though these results need to be verified with larger samples 
and more robust measurements. 
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APPENDIX A: OVERVIEW OF OUR STAR 
FORMATION LAW 

Here we give an overview of the star formation law used in 
our models. 

Following Blitz & Rosolowsky (2006) we assume that 
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Figure Al. Radial star formation rate surface density (SsFR (i?)) - gas surface density relations for 100 galaxies from our EFB model 
at rcdshifts 2 = (upper panels); 2 = 2 (middle panels) and 2; = 4 (lower panels). Each dot shows the star formation rate surface density 
and gas surface density in a single radial bin. Each galaxy has lO's of such radial bins. The three panels at each redshift correspond to the 
radial: atomic gas surface density (Satom(^J), left panels); molecular gas surface density (Snjoi(iJ), middle panels); and total gas surface 
density (Sgas(i?), left panels). The red lines show the theoretical asymptotic relations expected in the high pressure (Sgas — ^mol): ^i^d 
low pressure (Sgas ^ Satom) regimes. 
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star formation takes place in dense molecular gas, traced 
by the Hydrogen Cyanide molecule (HCN), with a constant 
star formation efficiency: 

SsFR eSF Sniol.HCN 



(Al) 



[Mo pc-2 Gyr-i] [Gyr"!] [Mq pc" 

where ?sf ~ 10 - ISGyr"^ (Gao & Solomon 2004, Wu 
et al. 2005). Expressing this equation in terms of the total 
gas content: 

SsFR esF Sg, 



[Mo pc-2 Gyr- 



[Gyr-i] [Mo pc" 



■/mol7?,HCN, (A2) 



where 72,hcn = Smoi.HCN/Smoi is the ratio between the 
dense molecular gas (as traced by HCN) and the total molec- 
ular gas, and /moi is the molecular gas fraction. 

The fraction of gas that is molecular, /moi, can be de- 
fined in terms of the mass ratio between molecular and 
atomic gas, T^-moi by 



(A3) 



"R-mol + 1 

Empirically Wong & Blitz (2002) and Blitz & Rosolowsky 
(2004; 2006) have argued that 7?.moi is determined to first 
order by the mid plane pressure, Pcxt. We adopt the relation 
from Blitz & Rosolowsky (2006): 

0.92±0.1 



4.3 ±0.6 X 10* 



(A4) 



where k is Boltzmann's constant, and Pcxt/k is in cgs units. 
For a gas plus stellar disk the mid plane pressure is given, 
to within 10%, by (Elmegreen 1993) 



-GEg 



V O" star / 



(A5) 



where cTgas and (Tstar are the velocity dispersions of the gas 
and stellar disk, respectively. For simplicity, we will assume 



0.1. 



Based on the arguments and references in Blitz & 
Rosolowsky (2006) we adopt the following fitting function 
for T^HCN 

0.4 



7^HCN = 0.1 1 -I- 



[200 Mo pc 



(A6) 



In the low pressure regime, /moi — 7?.moi oc Eg 



and T^HCN — 0.1, and thus in terms of total 



gas or atomic gas Eq. (|A1|) asymptotes to 



jSFR 



isF 1.37 X 10' 



[Mopc-2Gyr-i] [Gyr-i] [Mq pc-2] 

X 



[IMopc-^V' ■ ^^^^ 
Or in terms of the molecular gas, Eq. (|A1|I asymptotes to 

EsFR. _ £SF OA Emol / . „N 

[Mo pc-2 Gyr-1] " [Gyr-i] [Mq pc-2] [Mo pc'^] ' ^ ' 

In the high pressure regime, /moi — 1, Egas ^ Emoi, and 
TZhcn oc Emoj, and eQ. (|Al|) asymptotes to the familiar SK 
relation 

EsFR _ gsF 0.0120 

[Mo pc-2 Gyr-1] " [Gyr-i] [Mo pc-2] 

1.4 

(A9) 



[1 Mo pc- 



Furthermore, with esF ~ 13 Gyr^^, we recover the coefficient 
of esF — 0.16 Mo pc~^ Gyr~^ of the standard SK relatiorQ 
We implement the star formation recipe given by 
Eg. ljAlll as follows. At each time step and annulus in the 
disk, we calculate the star formation rate. Then we use the 
following approximation (valid for times steps small com- 
pared to the star formation time scale) to calculate the mass 
of newly formed stars 



AM,tar(ii) = A{R) EsFR(i?,t) At, 



(AlO) 



with A the area of the annulus, and At the time step interval. 
Note that At is not a constant in our model, it varies from 
~ 4 Myr at high redshifts to ~ 100 Myr at low redshifts 
(green dotted line in Fig. IB1[I . 

Fig. lAll shows the radially resolved star formation law 
from 100 randomly selected galaxies from our model with en- 
ergy driven outflows (EFB) at redshifts z = 0, 2, &4. Each 
dot corresponds to the star formation rate surface density 
and gas surface density of a single radial bin in a single 
galaxy. Each galaxy has lO's of such bins. The red lines 
show the asymptotic behaviour expected from our model 
in the high and low pressure regimes. In the high pres- 
sure regime Egas — Emoi, and so the star formation law 
(EsFR oc Egas) for the molecular and total gas have the 
same slope of n = 1.4 (dotted lines). In the low pressure 
regime the star formation law for the molecular gas has a 
slope of n = 1 (solid lines), while that of the atomic and 
total gas have a slope of n = 2.84 (dashed lines). The pan- 
els on the left show that there is a maximum surface density 
~ 15Mo pc~^ that atomic gas reaches (see Button (2009) for 
examples of atomic and molecular gas density profiles from 
our models). This causes the relation between Esfr and 
Eatom to steepen around Eatom — lOMo pc~^, and for the 
slope to reverse sign for Esfr ^ lOMopc"^. This regime 
corresponds to small galactic radii, where the atomic gas 
densities decrease towards the center of the galaxy, but the 
molecular gas densities and star formation rate densities are 
increasing. The features of our star formation laws closely 
resemble the observed relations between local star forma- 
tion rate surface density and atomic, molecular, and total 
gas density as measured by Bigiel et al. (2008), as well as 
the theoretical model of Krumholz, McKee, & Tumlinson 
(2009). 



APPENDIX B: HOT MODE VS. COLD MODE 
ACCRETION 

It has recently become clear that the way we model the 
accretion of gas into galaxies (via a cooling fiow of shock 
heated gas), is not entirely correct. Birnboim & Dekel (2003) 
and Dekel & Birnboim (2006) argued that in haloes below a 
characteristic mass galaxies accrete most 

of their gas via "cold fiows" , in which the gas is never heated 
to the virial temperature of the halo. This has subsequently 
been confirmed in cosmological, hydrodynamical simulations 



^ Note that the coefficient of egF = 0.25 Mo pc"'^ Gyr~^ from 
Kennicutt (1998) does not account for Helium in the gas masses. 
Correcting for Helium with a factor of 1.36 results in egF = 
0.16 Mq pc-2 Gyr-i. 
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Figure Bl. Several relevant timescales in our model as function 
of redshift: Hubble time (blue short-dashed); free fall time (red 
solid); cooling times for three halo masses (black long-dashed); 
and numerical time step (green dotted). Gas accretes onto the 
central galaxy on a timescale that is the maximum of the cooling 
time and the free fall time. At high redshifts the free fall time is 
the limiting timescale, and thus there is a delay of a few 100 Myr 
between when gas enters the halo and when it accretes onto the 
central galaxy. For massive haloes at low redshifts, the cooling 
time exceeds the free-fall time, and may even exceed the Hubble 
time. 

(e.g., Keres et al. 2005). According to Dave (2008), the defin- 
ing features of cold mode accretion are that it is (i) rapid, (ii) 
smooth, and (iii) steady. Rapid because it is limited by the 
free fall time, rather than the cooling time, smooth because 
most (more than 90%, Dekel et al. 2009) of the accretion 
occurs in small lumps and not through major mergers, and 
steady because it is governed by the gravitational potential 
of the slowly growing halo. 

What differentiates the cold mode accretion model from 
the classical hot mode accretion model? In the hot mode 
model, gas that enters the halo is shock heated to the virial 
temperature. The hot gas then cools radiatively, and the 
cold gas accretes onto the central galaxy on a timescale that 
is the maximum of the cooling time or the free fall time. In 
cold mode accretion there is either no virial shock, or cold 
streams penetrate the hot halo, and gas accretes onto the 
central galaxy on a free fall time. 

Fig. lBll shows the cooling times (tcooi, black long-dashed 
lines) and free fall times (iff, red solid line) versus redshift 
for three halo masses (Mvir = lO^^ ^'^ ^'^^ ^/i-^A/o) in our 
adopted cosmology. This shows that for a given halo, both 
the cooling time and free fall time decrease to higher red- 
shift, but that the cooling time decreases faster. Thus at 
high redshifts, the free fall time is the limiting factor in the 
cold gas accretion rate. For typical halo masses of star form- 
ing galaxies (Mvir — 10^^~^'^Mq) the cooling time is shorter 
than the free fall time for most redshifts. 

The impact of this competition between the cooling 
time and free fall time on the accretion rates is shown in 



Fig. IB2I This shows the accretion ratio vs. halo mass at 
redshifts z = to z = 5. The accretion ratio is defined as 
the ratio between the accretion rate of cold gas onto the 
galaxy, Mcooi, and the baryon accretion rate into the halo, 
Mbar, at that redshift. As expected from Fig. IBll at low 
halo masses and high redshifts haloes are in the fast cooling 
regime (tcooi < iff), and thus the accretion ratio is close to 
unity (deviations from unity arise because of the exponen- 
tial mass accretion histories adopted here). Hence, although 
we do not explicitly model cold accretion flows, effectively 
the model yields cold gas accretion rates that are indistin- 
guishable from those expected from cold flows. 

Above a halo mass of Mvir,z — 5 x IO'^^Mq, the haloes 
are in the slow cooling regime (tcooi > iff), and the accre- 
tion ratio drops sharply. It has been argued that in haloes 
above this mass range at high redshifts (2 > 2), the cold gas 
accretion rates are dominated by cold streams penetrating 
the hot halo (Dekel et al. 2009) . This aspect of gas accretion 
is not captured by our cooling/accretion treatment, which 
may impact the formation of massive galaxies at z > 2. 
Apart from this specific case, the classical hot mode cooling 
model that we adopt in this paper, reproduces the defining 
features of cold mode accretion: rapid, smooth, and steady 
(e.g. Dave 2008). We note that while we refer to cold accre- 
tion as smooth, it can also be described as clumpy in the 
sense that the accretion is made up of many small lumps, 
rather than a diffuse structureless medium. 



APPENDIX C: WHY ARE HIGH REDSHIFT 
DISKS DENSER? 

As discussed in H3.3.3I the evolution of the densities of gas 
disks in our model is stronger than predicted by the simple 
halo density scaling. To understand what causes this devia- 
tion we first discuss what causes the evolution in the surface 
densities of the baryonic disks in our model. Since we assume 
that the disks are in centrifugal equilibrium, the evolution 
in the density of the disk is determined by the evolution of: 
1) the concentration of the halo, c; 2) the ratio between the 
disk mass and total mass, rUgai; and 3) the ratio between 
the spin parameter of the disk and the spin parameter of 
the halo (i.e. the ratio between the specific angular momen- 
tum of the disk and the specific angular momentum of the 
halo), jgai/mgai. Higher c, higher mgai and lower jgai/rrigai 
all result in higher density disks. 

To isolate the effects of evolution in c from that of mgai 
and jgai/?Wgai, we introduce a new model, labeled NFC, in 
which we have turned the feedback and cooling calculations 
off. In this model we assume that the ratio between the 
baryonic mass and virial mass, mgai = 0.04, independent of 
redshift. To implement this, at each time step, we set the 
baryonic mass that accretes onto the galaxy to be rrigai times 
the total mass (baryons plus dark matter) accreted into the 
virial radius during that time step. In this model we also 
assume that the spin parameter of the halo is independent 
of redshift, with a median value of A = 0.035, and that the 
spin parameter of the disk is equal to the spin parameter of 
the halo, i.e. jgai/mgai = 1. 

The evolution of the surface density of the baryonic 
disks in the NFC model is shown as open circles in the lower 
left panels of Fig. ICll This evolution follows the virial scal- 
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Figure B2. Ratio between the accretion rate of cold gas onto the central galaxy, Meool and the baryon accretion rate into the virial 
radius, M^ar (the accretion ratio) as a function of the virial mass at redshift zero (left-hand panel) and as a function of the virial mass 
at redshift z (right-hand panel). The halo mass of Afvir,z — 5 X 10^^ Mq marks the transition between fast cooling (tcool < iff) a^nd slow 
cooling (tcool > iff )- Hence, our hot- mode cooling model yields cold gas accretion rates that are indistinguishable from those expected 
from cold flows. 



ing for redshifts z > 1, but for redshifts ^ < 1 the evolu- 
tion is weaker. The weak evolution at low redshifts is due 
to the evolution in the halo concentration (see Somerville 
et al. 2008). 

The middle and right panels in Fig. ICll show the evo- 
lution of jgai/'Tigai and JTigai. By construction the NFC 
model (open circles) shows no evolution in these parame- 
ters. The other three models result in jgai/?Tigai decreasing 
with increasing redshift. For the NFB model the evolution 
in jgai/?7igai is due to the inside-out nature of the cooling 
and accretion of the hot halo gas in our model. As shown in 
Fig. IBll at high redshifts the cooling time of the hot halo 
gas is less than the free fall time, and thus there is a delay 
between when baryons enter the halo, and when they reach 
the central galaxy. At z ~ 6 roughly 60% of the baryons that 
accreted into the halo have cooled and reached the central 
galaxy. Since the mass of the halo is more centrally concen- 
trated than the angular momentum, by only accreting the 
inner 60% of the hot halo baryons, the gas that reaches the 
disk has lower specific angular momentum than the halo. By 
redshift z = 0, rrigai — /bar, and jgai/?Tigai — 1- For the NFB 
model jgai/rrigai decreases to higher redshifts, which should 
result in higher density disks. However, nigai also decreases 
to higher redshifts, which counteracts this effect. This ex- 
plains why the evolution of the disk density in the NFB 
model follows the evolution of the virial scaling. 

At high redshifts all three models with cooling (EFB, 
MFB, NFB) result in similar jgai/?7tgai and rrigai. This is be- 
cause outflows are inefficient in high redshift galaxies due to 
the deep potential wells. However, outflows do occur at high 
redshifts, which result in the EFB model having lower mgai 
and higher jgai/»7igai than the MFB and NFB models. Out- 
flows result in higher jgai /rrigai because gas is preferentially 
removed from small radii (where the stars are preferentially 
forming due to the density dependent star formation law). 



At lower redshifts, the galaxies are less dense, escape veloc- 
ities are lower, and thus outflows are more efficient. For the 
EFB model this results in rrigai decreasing with decreasing 
redshift. Thus the combined effects of lower specific angu- 
lar momentum and higher mass fractions in higher redshift 
disks result in the disk densities increasing faster than the 
virial scaling. 

Now we return to the evolution of the gas densities. The 
evolution of the gas density for the NFC model is shown in 
the upper left panel in Fig. ICll This shows a similar, but 
slightly stronger evolution than the virial scaling. This de- 
viation is qualitatively explained by the spin parameter of 
the gas decreasing with increasing redshift and the gas mass 
fraction increasing with increasing redshift (upper middle 
panel). The evolution of the gas densities in the NFB, MFB 
and EFB models is stronger than that of the NFC model 
because the specific angular momentum of the disk evolves 
faster with redshift in these models. Thus in our model the 
gas densities evolve more strongly than the virial scaling be- 
cause the spin parameters of the gas (and baryons) decrease 
to higher redshifts. 



© 2009 RAS, MNRAS 000.[Tlf23l 



24 Button et al. 




Iog,o(1+z) log,o(1+z) log,o(1+z) 

Figure CI. Upper panels: Evolution in the effective surface density of the cold gas (left); evolution in the ratio between the specific 
angular momentum of the cold gas to that of the halo (middle); evolution in the cold gas mass fraction (right). Lower panels: same as 
upper panels but for cold gas+stars. The deviation from the virial scaling (red long dashed lines) is driven by evolution in the ratio 
between the specific angular momentum of the disk and the specific angular momentum of the halo. See text for further details. 
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